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Abstract Non-thermal components are key ingredients for understanding clusters of 
galaxies. In the hierarchical model of structure formation, shocks and large-scale tur- 
bulence are unavoidable in the cluster formation processes. Understanding the ampli- 
fication and evolution of the magnetic field in galaxy clusters is necessary for mod- 
elling both the heat transport and the dissipative processes in the hot intra-cluster 
plasma. The acceleration, transport and interactions of non-thermal energetic parti- 
cles are essential for modelling the observed emissions. Therefore, the inclusion of the 
non-thermal components will be mandatory for simulating accurately the global dy- 
namical processes in clusters. In this review, we summarise the results obtained with 
the simulations of the formation of galaxy clusters which address the issues of shocks, 
magnetic field, cosmic ray particles and turbulence. 

Keywords cosmology; theory, large-scale structure of universe; acceleration of 
particles; hydrodynamics; magnetic fields; method; numerical, N-body simulations 



1 Introduction 

Within the intra-cluster medium (ICM), there are several processes which involve non- 
thermal components, such as the magnetic field and cosmic rays (CRs). At first glance, 
many of these processes can be studied numerically in simplified configurations where 
one can learn how the individual processes work in detail. For example, there are 
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investigations which co nsider maKneto-hyd ro-dynamical (MHD) simulations of cloud- 
wind interactions (see iGreeori et al. 2000l. and referen c es therein) or simula te the 
rise of relic radio bubbles (see Jones fc De Youn"3 boosi : [Reynolds et alj boosl . and 
references herein). Such work usually focuses on the relevance of CRs, turbulence and 
the local magnetic fields within and around these bubbles. In this review, we will instead 
concentrate on simulations of non-thermal phenomena, which aim at understanding 
the relevance of these phenomena for galaxy cluster properties or at unveiling possible 
origins of the non-thermal radiation. So far, firm evidence for the presence of non- 
thermal emission (at radio wavelengths) and for the presence of cosmologically relevant 
extended magnetic fields has been found only in galaxy clusters; in filaments only weak 
limits have bee n found so far. For a m ore detailed discussion see th e reviews on radio 
observations bv lRephaeli et ah! 1 20081 ) - Chapter 5, this volume, and lFerrari et allboosi 
- Chapter 6, this volume. 



2 Possible origins of non-thermal components within the large-scale 
structure (LSS) 

The origin of magnetic fields and of CRs in galaxy clusters is still under debate. The 
variety of possible contributors to magnetic fields ranges from primordial fields, battery 
and dynamo fields to all the classes of astrophysical objects which can contribute with 
their ejecta. The latter possibility is supported by the observatio n of metal enrichmen t 
of the ICM, which is due to the ejecta of stars in gal axies, see IWerner et al.l boOSi - 
Chapter 16, this volume, and lSchindler fc DiaferigbOOSl - Chapter 17, this volume. The 
magnetic fields produced by all these contributors will be compressed and amplified by 
the process of structure formation. The exact amount of amplification and the resulting 
filling factor of the magnetic field will depend on where and when the contributor is 
thought to be more efficient. 

CRs are produced at shocks, both in supernova remnants and in cosmological ac- 
cretion and merger shocks. Furthermore, active galactic nuclei (AGN) and radio galax- 
ies contribute to the CR population. CR protons can also produce CR electrons via 
hadronic reactions with the thermal plasma, leading to the so-called secondary CRs. 



2.1 Possible origins of magnetic fields 

Magnetic fields can be produced either at relatively low re dshift (z ~ 2 — 3) or at 
high redshift {z > 4). In the former case, galactic winds (e.g. IVolk fc AtovanlbOOCj ) or 
AGN ejecta (e.g. lEnfilin et al.|[l997l : iFurlanetto fc LoeblboOll and references therein) 
produce magnetic fields 'locally', e.g. within the proto-cluster region. In the latter case, 
the magnetic field seeds can also be produced by an early population of dwarf starburst 
galaxies or AGN before galaxy clusters form gravitationally bound systems. 

One of the main arguments in favour of the 'low-redshift' models is that the high 
metallicity observed in the ICM suggests that a significant enrichment occurred in the 
past due to galactic winds or AGN. Winds and jets should carry magnetic fields to- 
gether with the processed matter. It has been shown that winds from ordinary galaxies 
give rise to magnetic fields which are far weaker than those observed in galaxy clus- 
ters, whereas magnetic fields produced by the ejecta of starburst galaxies can be as 
large as 0.1 fiG. Clearly, this class of models predicts that magnetic fields are mainly 
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concentrated in and around galaxies and within galaxy clusters. Note that if the mag- 
netic pollution happens early enough (around z ~ 3), these fields will be amplified not 
only by the adiabatic compression of the proto-cluster region, but also by shear flows, 
turbulent motions, and merging events during the formation of the cluster. Shocks 
are expected to be produced copiously during the non-linear stage of the LSS forma- 
tion (see below). Recent detailed studies of shock propagation re vealed the presence 
of specific instabilities driven by energetic accelerated particles feell fc Lucek' ^OOll : 
IVladimirov et aklbood ). Such instabilities result in a strong, non-adiabatic amplifica- 
tion of an upstream magnetic field seed which converts an appreciable fraction of the 
shock ram pressure into magnetic fields. 

In the 'high-redshift' models of the magnetic field generation, the strength of the 
field seed is expected to be considerably smaller than in the previous scenario, but 
the adiabatic compression of the gas and the shear fiows driven by the accretion of 
structures can give rise to a considerable amplification of the magnetic fields. Several 
mechanisms have been proposed to explain the origin of magnetic field seeds at high 
redshift. Some of them are similar to those discussed above, differing only in the time 
when the magnetic pollution is assumed to take place. In the present class of models the 
magnetic field seeds are supposed to be expelled by an early population of AGN or dwarf 
starburst galaxies at a redshift between 4 and 6 ( Kronberg et al.lll999l ): this process 
would magnetise a large fraction of the volume. Recently, the vali dity of such a scenari o 
has been confirmed by a semi-analytic modelling of galactic winds i Bertone et al.|[2006l ). 

Alternative models invoke processes that took place in the early universe. Indeed, 
the presence of magnetic fields in almost all the regions of the universe suggests that 
they may have a cosmological origin. In general, most of the 'high-z models' predict 
magnetic field seeds filling the entire volume of the universe. However, the assumed 
coherence length of the field crucially depends on the details of the model. While 
scenarios based on phase transitions give rise to coherence lengths which are so small 
that the corresponding fields have probably been dissipated, magnetic fields generated 
at neutrino or photon decoupling are thought to have a much higher chance of surviving 
until the present time. Another (speculative) possibility is that the field seed was 
produced during i nfiation. In this case, t he co herence length can be as large as the 
Hubble radius. See lCrasso fc RubinsteinI ( 200 ih for a review. 

Magnetic field seeds can also be p roduced by the so-called Biermann battery effect 
( Kulsrud et Zl 1 19971 : IRvu et al.l[l998l ). The idea here is that merger/accretion shocks 
related to the hierarchical structure formation process give rise to small thermionic 
electric currents which, in turn, may generate magnetic fields. The battery process has 
the attractive feature of being independent of unknown physics at high redshift. Its 
drawback is that, due to the large conductivity of the ICM, it can give rise to very tiny 
magnetic fields, of order 10~^^ G at most. One therefore needs to invoke a subsequent 
turbulent dynamo to boost the field strength to the observed level. 

Such a turbulent amplification, however, cannot be simulated numerically yet, mak- 
ing it quite difficult to predict how it would proceed in a realistic environment. It is 
clear that one expects the level of turbulence to be strongly dependent on the environ- 
ment, and that it should appear mostly in high-density regions like collapsed objects. 
While energetic events, such as mergers of galaxy clusters, can easily be effective at 
driving the required levels of turbulence, it is harder to understand how turbulent am- 
plification would work in relatively quiet regions like filaments. Here, magnetic fields 
can be rather amplified by the shocks which originate during the LSS formation, as 
mentioned at the beginning of this section. Lacking a theoretical understanding of the 
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turbulent amplification, it is therefore not straightforward to relate the very weak field 
seeds produced by the battery process to the magnetic fields observed today. Attempts 
to construct such models, based on combining numerical and analytical computations, 
have not yet been reported to successfully reproduce the observed scaling relations of 
the magnetic fields in galaxy clusters. 



2.2 Possible origins of CRs 

There are many possible contributors to the CR population in the ICM. AGN and 
stellar activity in cluster galaxies are believed to be a significant source of CRs in the 
ICM. Furthermore, cluster accretion shocks of high Mach numbers (typically above 
100) and of sizes of megaparsec or larger are wid ely accepted to be possible CR sources 
with nuclei accelerated up to ~ 10^ GeV (e.g. iNorman et al.| [l995). Cluster mergers 
generate internal shocks (of moderate Mach numbers, ty pically below 4) which provide 
most of the ICM gas heating (see e.g. iKang et al.ll2005[ ). but are also likely to convert 
a non-negligible fraction (> 10 %) of their power into CRs. 

AGN outfiows dissipate their kinetic energy and Poynting fiuxes into the ICM 
providing additional non-gravitational I CM heating and, thus, a plausible feedbac k 
solution to the cooling flow problem (e.g. lBinney fc Taborll 19951 : IChurazov et al.ll2003 ) 



Recent XMM-Newton observations of the cluster MS 0735-1-7421 (jCitti et all boOTi ') 
revealed that a powerful AGN outburst deposited about 6x10®^ erg into the ICM 
outside the cooling region of radius ~ 100 kpc. Relativistic outfiows of AGN are likely to 
be an essential source of super-thermal particles in clusters. A powerful relativistic AGN 
jet could deposit up to 10^^ erg into a relativistic particle pool during a duty cycle of 
about ~ 50 Myr. A potential source of energetic particles are sub-relativistic ions with 
energy 100 MeV that a re evaporated fro m an AGN, which accretes mass in the ion- 
supported tori regime fe.g. lRees et al.|[l982h . A certain fraction of the energetic particles 
could escape the fiows and avoid fast cooling due to collective effects in the central parts 
of the cluster. Re-acceleration of that non-thermal population by inner shocks inside 
the cluster can provide a long-lived non-thermal component that contributes to the 
total pressure of the ICM. 

S tar formation act ivity in galaxies is another source of CRs in clusters of galaxies 
(e.g. IVolk et al.l Il99d . and references therein). The combined action of supernovae 
and winds of early-type stars leads to the formation of a hot. X-ray emitting, slowly 
expanding superbubble filled with large-scale (tens of parsecs) compressible MHD mo- 
tions (e.g. rarefactions) and shocks. When the starburst event is energetic enough, the 
superbubble may expand beyond the disk of the parent galaxy and produce a super- 
wind that supplies the intergalactic medium with metals ejected by supernovae and 
CRs. The superwinds have been seen with Chandra in some starburst galaxies: M 82, 
Arp 220 and NGC 253. In the superwind scenario, the CR ejection is tightly connected 
with the metal enrichment of the cluster (|Schindler fc Diaferiolliooj - Chapter 17, this 
vol ume). 

iBvkovl l|200ll ) argues that non-thermal particle acceleration can favourably occur 



in correlated supernovae events and powerful stellar winds with great energy release, 
which generate interacting shock waves within superbubbles. The acceleration mech- 
anism provides efficient creation of a non-thermal nuclei population with a hard low- 
energy spectrum; this population can transport a substantial fraction (~ 30%) of the 
kinetic energy released by the supernovae and by the winds of young massive stars. 
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Fig. 1 The simulated efficiency of conversion of shock and bulk MHD motions ram pressure 
power to the CRs at different stages of cluster evolution. The efficiency is defined as a volume 
averaged fraction of the cluster large-scale motions power transferred to CRs. The simulation 
was done for a model of CR acceleration by an ensemble of shocks and MHD-turbulent motions 
described in^Bykov ([2001), adopted for clusters of galaxies. The MHD turbulence in a cluster 
of a Mpc size was initiated with the Gaussian bulk velocity dispersion of 2 000 km s~^ at a 
scale of 100 kpc. 



A bright phase in the galaxy evolution can be the source of the rehc CRs in clusters. 
Energetic nuclei can b e stored in cluster magnetic fields for several Hubble times (e.g. 
iBerezinskv et al.lll997l '). The presence of these nuclei produces a diffuse fiux of high- 
energy 7 and neutrino radiation du e to the interaction of the CRs with the ICM 
( Volk et al.lll996l : lEnfilin et al.lll997l l. The resuhing flux depends on the ICM baryon 
density. 

A physical model of particle acceleration by the ensemble of inner shocks in the 
ICM could be similar to the superbubble model. The sub-cluster merging processes 
and the supersonic motions of dark matter halos in the ICM are accompanied by 
the formation of shocks, large-scale flows and broad spectra of MHD-fluctuations in 
a tenuous intra-cluster plasma with frozen-in magnetic fields. Vortex electric fields 
generated by the large-scale motions of highly conductive plasma harbouring shocks 
result in a non-thermal distribution of the charged nuclei. The free energy available for 
the acceleration of energetic particles is in the ram pressure of the shocks and in the 
large-scale motions. 

The most studied way to transfer the power of the MHD motions to the energetic 

parti cle p opulation is the Fe rmi-type acceleration (see e.g. the review bv lBlandford &: Eichled 
Il987l andl&yk ov et al.ll2008l - Chapter 7, this volume). An important ingredient of the 
energetic particle acceleration by shocks and large-scale MHD motions is the presence 
of small-scale MHD turbulence, which is necessary to scatter relativistic particles and 
to make their pitch- angle isotropic. The scale of the fluctuations required for the reso- 
nant scattering of a particle of energy ^ 1 GeV is about 3x 10^^B"i cm, where B-r 
is the local mean magnetic field in pG. The scale is some 10-11 orders of magnitude 
smaller than the basic energy scale of the system; thus the origin and the maintenance 
mechanisms of such small scale turbulence is a serious problem. 
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Fig. 2 The temporal evolution of the particle distribution function in a cluster of a Mpc size 
and a Gaussian bulk velocity dispersion of about 2 000 km s~^ at a scale of 100 kpc. The 
curves labelled as 1, 2, 3 correspond to 0.1, 0.3 and 1 Gyr age of the cluster. Mono-energetic 
proton injection at energy 10 keV (po is the injection momentum). The injection efficiency 
(see Fig.[TJ was assumed to be about 10""^. The power conversion efficiency for the model was 
already illustrated in Fig. [T] 



In non-linear models of particle acceleration by strong MHD shocks the presence 
of turbulence could be supported by t he CR instabilities them selves (see e.g. recent 
non-linear Monte-Carlo simulations by Iviadimirov et al.ll200^ ). Direct MHD cascade 
of energy to small scale turbulence is also p ossible, but the small scale fluctuations 
are highly anisotropic in that case (see e.g. lBiskamTjl2003h . Moreover, the cascade 
properties are still poorly known at the scales close to the Coulomb mean free path 
where the transition from the colhsional to the collisionless plasma regime occurs. 



The particle distribution within a system with multiple shocks and with large-scale 
plasma motions is highly intermittent, because it has strong peaks at shocks and a non- 
steady mean CR background. By using the kinetic equations given by |Bvkovll|200ll ). one 
can construct a model of the temporal evolution of the particle distribution function, 
which accounts for the non-linear effect of the reaction of the accelerated particles on 
the shock turbulence inside the cluster. In Fig. [T] we show the efficiency of a conversion 
into CRs of the power of shocks and the bulk plasma motions of a scale of about 
100 kpc. One may note that, while the efficiency could be as high as 30 % for some 
relatively short period of time, the efficiency is typically ~ 10 % for most of the time. 
Fig. [2] shows the temporal behaviour of the particle spectra on a Gyr time scale. 



Non-thermal particles and magnetic fields could contribute to the total pressure of 
the ICM. In the particular model described above a substantial energy density could 
be stored in protons of energy below 200 MeV. Protons of that energy level produce 
very little high energy emission and therefore it is rather difficult to constrain that 
component directly from the observational data which are currently available. 
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3 Shocks in the LSS and clusters 

As described in tlie previous section, shocks can play a crucial role as source of non- 
thermal components in galaxy clusters. Therefore, a significant effort is spent to inves- 
tigate the distribution and the evolution of shocks during structure formation. 

In early simulations of the formation of cosmological str uctures, it was alr eady 
noticed that shocks are the main source of the ICM heating fsee lOuilis et al.lll998l . and 
references therein). In numerical simulations, two classes of shocks can be identified. 

The first one is the accretion shocks surrounding the forming objects, namely fila- 
ments or galaxy clusters. In the case of galaxy clusters, most shocks are quasi-spherical 
and their position with respect to the centre of the cluster grows with time (shock trav- 
elling outwards). These shocks can be associated with the accretion of diffuse matter 
which is relatively cold. These shocks provide a first step towards the virialisation of 
an initially supersonic accreting flow. The radius of the accretion shock (for a quasi- 
spherical cluster at z = 0) is typically two times the virial radius. As the typical 
temperature and density of the upstream gas is quite low, such shocks can have very 
large Mach numbers (S> 10). 

The second class of shocks in galaxy clusters arises from the merging of substruc- 
tures. The gas in gravitationally bound substructures is much denser than the diffuse 
ICM; therefore it cannot be stopped in the accretion shock and it continues to move 
with the substructure through the gas, which is mostly virialised. The substructures, 
which have velocities larger than the sound velocity of the hot ICM, drive the so-called 
merger shocks. These shocks are initiated within the ICM which has been already 
shock-heated, and therefore they typically have modest Mach numbers (below !«4). 
The merger shocks propagate within the dense ICM and are likely to be the main 
contributor to the global energy dissipation within the ICM (see lPfrommer et ablbooil . 
and references therein). As these shocks are primarily driven by gravitational attrac- 
tion, non- gravitational processes (like cooling or feedback by s tar formation proc esses) 
do not significantly affect the global energy dissipation (see iKang et al.ll2007l ). The 
ram pressure dissipation process in shocks within the coUisionless plasma is non-trivial 
because of the role of long-lived non-thermal co mponents, namely magnetic fields, 
MHD-waves and energetic particles (see for details fevkov et ak I boOSi - Chapter 7, this 
volume) . The complex multi-scale structure of a shock modified by energetic CRs (e.g. 
where the back reaction of the CRs on the thermal plasma is strong enough to change 
the thermodynamical states of the plasma in addition to the shock itself) cannot be 
resolved properly in LSS simulations. 

The efficiency of energising CRs by shocks strongly depends on a number of factors. 
Among the most important factors, we have (i) the upstream magnetic field of the 
shock, (ii) the plasma parameter (3 = M^/Ms (ratio of thermal to magnetic pressure, 
here expressed by the Mach numbers Aia and A^s with respect to the Alfven and sound 
velocity respectively), (iii) the inclination of the upstream field of the shock, and (iv) 
the presence of a pre-existing CR population produced by the accretion shock or in 
star-forming regions. 

The distribution of the Mach numbers of the shoc ks within the cosmological struc- 



ture can be derived from semi-analytical modelling (see lCabici fc Blasilbood : IPavhdou fc Fields! 



20061. and references therein), as well as d irectly fro m cosmological simulations (see 



Miniati et al. 2000: Ryu et al. 2003; Pfromn ^eret al.i r2006: Kang ct al. 2007', and ref- 
erences therein). An example is shown in Fig. [3] There are significant differences in 
the strength and the distribution of the shocks found in different works: this fact has 
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Fig. 3 Shocks detected in a cosmological simulation of a cubic volume with a side length of 
100 Mpc, resolved with 512^ grid points. The figure shows the kinetic energy flux per unit 
comoving volume through surfaces of external and internal shocks with Mach numb er between 
logM and logM+d(logM) at different redshifts (solid, z = to dotted, z = 2). From lRvu et al.l 
l|2003D . 



strong consequences on the predicted amount and the relative importance of CRs and 
their energy with respect to the thermal energy of the ICM. Note that part of the 
confusion, specially in the early work, is caused by a lack of resolution with in the sim- 
ulations, which prevents int ernal shocks from being properly resolved (see IRyu et ahl 
booi IPfrommer et aLlliooH ): further confusion originates from the influence of other, 
non-therma l processes (like preheating) on the infe rred Mach number of the external 
shocks (see IPfrommer et aLlbood : iKang et al.ll2007l l. Moreover, differences in the un- 
derlying hydrodynamic simulations as well as in the shock detection, mostly done in 
a post processing fashion, can contribute to differences in the results. Nevertheless, 
there is a qualitativ e agreemeii t in th e ca se of simulated internal shocks between the 
simulations done bv lRvu et ai] (|2003l ) and IPfrommer et al.l (|2006l ). 



4 Evolution of magnetic fields in simulations 

4.1 Local amplification of magnetic fields 

A very basic process which amplifies magnetic fields is related to the Kelvin-Helmholtz 
(KH) instabilities driven by shear flows, which are common during the formation of 
cosmic structures. iBirk et all (|l999l l performed a detailed study of such an amplification 



within the outfiows of starburst galaxies, where the KH timescale should be « 4 x 
10^ years (see Fig. By using a Cartesian resistive MHD code they found that 
the amplification factor of the magnetic field mainly depends on the initial ratio of the 
magnetic to the kinetic energy and only mildly depends on the assumed resistivity. They 
concluded that such a process could indeed explain why the magnetic field observed 
in the halo of starburst galaxies is significantly higher than what is expected from 
the magnetic fields observed within galactic disks. When applied to a cluster core 
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Fig. 4 Resistive, MHD simulation of a flow which is KH instable at the interface separating 
the left and right halve of the simulated volume. Shown is the magnetic field (y-component) 
in units of the initial magnetic field after evolving the KH instability (inlay) and the time 
evolution of the (relative) magnetic field amplitude for different values of the chosen resistivity 
and initial magnetic field. Time is in units of the Alfvenic transit time (4000 year). Triangles 
and stars are for an initital magnetic field of 10 ^G, crosses and diamonds for a ten times larger 
field. Stars and diamonds have a magnetic Reynolds number of 100, tr iangles and crosses have 
a ten times larger magnetic Reynolds number. From ISirk et all ||1999D . 



environment, the KH timescale turns out to be 10^ years; this timescale makes the KH 
instability an interesting process to amphfy weak magnetic fields. 

It is quite expected that merger events and accretion of material onto galaxy clus- 
ters will drive significant shear flows within the ICM. Extensive MHD simulations 
of single merging events performed using the Eulerian code ZEUS I Stone &: NormanI 
Il992a b ) demonstrated that merger events are effective at amplifying magnetic fields 
(jRoettiger et aLlll999bl ). In particular these authors found that the field initially be- 
comes quite filamentary, as a result of stretching and compression caused by shocks 
and bulk flows during infall; at this stage only a minimal amplification occurs. When 
the bulk flow is replaced by turbulent motions (e.g., eddies), the field amplification is 
more rapid, particularly in localized regions, see Fig.O The total magnetic field energy 
is found to increase by nearly a factor of three with respect to a non-merging cluster. 
In localised regions (associated with high vorticity), the magnetic energy can increase 
by a factor of 20 or more. 

A power spectrum analysis of the magnetic energy showed that the amplification is 
largely confined to scales comparable to or smaller than the cluster cores: this indicates 
that the core dimensions define the injection scale. It is worth taking notice that the 
previous results can be considered as a lower limit on the total amplification, because 
of the lack of resolution. Furthermore, it is quite likely that a galaxy cluster undergoes 
more than one of these events during its formation process, and that the accretion of 
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Fig. 5 The evolution of (the logarithm of) the gas density (left column), the gas tempera- 
ture (central column), and (the logarithm of) the magnetic pressure (right column) in two- 
dimensional slices taken through the core of a cluster in a major merger phase in the plane of 
the merger. Each row refers to different epochs: t = (i.e. the time of the core coincidence), 
t = 1.3, t = 3.4, and t = 5.0 Gyr, from top to bottom. Each panel is 3.75 X 3.75 Mpc. The 
merging subcluster enters from the right. From iRoettiger et al.l dlQQQbT) . 



smaller haloes also injects turbulent motions into the ICM: consequently, the magnetic 
field amplification within galaxy clusters will be even larger. A detailed discussion of the 
amplification of magnetic field in the cluster environment, using various simulations of 
driven turbulence, can be found in Subramanian ct al. (2006); these authors show that 
turbulent processes can provide reasonable strength and length scales of the magnetic 
fields in galaxy clusters. 
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Fig. 6 A two-dimensional cut through a cosmological box simulated bv lMiniati et al" 1 1I2OOII '). 
following the evolution of battery fields. The logarithm of the up-scaled magnetic field strength 
in Gauss at z = is shown. Taken from iSigl et al.. (.2004') . 



4.2 Magnetic fields from cosmological shocks 



Cosmological shocks, mainly the accretion shocks on cosmological objects like galaxy 
clusters and filaments, occur more frequently than the ones produced by individual 
merger events. They also can produce magnetic fiel ds through the so-called Biermann 
battery effect ( Kulsrud et al.|[l997l : IRyu et al.lll998l ). on which a subsequent turbulent 



dynamo may operate (see Sect. [5}. 

In such a scenario, the magnetic field is strongly correlated with the large-scale 
structure (see Fig. [B]). This means that the magnetic field within the filamentary 
structure could be even slightly higher than its equipartition value without violat- 
ing th£_Jwea^) U2£er limits of the rotation measure (RM) of quasars, as pointed out 



by IRvu et al.l 1 19981 ). It is worth pointing out that the arguments for the turbulent 



dynamo action, which could am plify the battery field seeds up to the /xG level, as pre- 
sented in lKulsrud et al.1 (1993), refer explicitly to regions which are close to collapse 



into galaxies. It has still to be proven that such arguments hold within proto-clusters or 
even cosmological structures, like sheets and filaments. In general, the time evolution 
of the magnetic field , as predicted by these simulations, strongly flattens around z ~ 3 



( Kulsrud et al]|l997l ): ;he evolution then leads to a magnetic field intensity which is 



relatively uniform on scales of tens of Mpc within the LSS around galaxy clusters (see 
Fig- El). Note that, so far, the synthetic rotation measures due to the magnetic fields 
predicted by up-scaling the battery fields have never been compared with the rotation 
measure observed on the scale of galaxy clusters. This might be partially motivated by 
the lack of resolution in the mock clusters within such non-adaptive, Eulerian simula- 
tions. 
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Fig. 7 The strength of the magnetic field as a function of baryonic overdensity within a cos- 
mological simulation. The straight Une shows the expectation for a purely adiabatic evolution, 
the s olid line gives the m ean field strength at a given overdensity within a cosmological simula- 
tion l lDolag et al.ll2005al ). While the latter is close to the adiabatic value in underdense regions, 
in clusters there is a significant inductive amplification due to shear fiows and turbulence; this 
amplification however saturates in the cluster cores. At any given density, a large fraction of 
particles remains close to the adiabatic expectation, as shown by the dotted line, which gives 
the median of the distribution at each density. 



4.3 Cosmological MHD simulations 



By using the GrapeMSPH code ( Dolag et alJ 1 19991 ) and assuming that a small ini- 



tial magnetic field seed exists before structure formation (see Sect. [Sjl, the first self- 
consistent simulations which follow the magnetic field amplification during the for- 
mation of galaxy cluster s within a cosmological environment have been performed by 
iDolag et al. I l|l999l , l2002l ). These simulations demonstrated that the contribution to the 



amplification of magnetic fields by shear flows (and by its induced turbulence) is signif- 
icant (see Fig. [7| ■ Therefore, for the flrst time, a consistent picture of the magnetic field 
in galaxy clusters could be constructed: the amplification predicted by the simulations 
was capable of linking the predicted strength of the magnetic field seed (see Sect. [2] 
and references therein) at high redshift (z « 3 and higher) to the observed magnetic 
field strength in galaxy clusters in the local universe. 

Furthermore, the simulations predicted that the final structure of the magnetic 
field in galaxy clusters refiects the process of structure formation, and no memory on 
the initial magnetic field configuration survives: this relaxes the constraints on models 
for magnetic field seeds. In general, such models predict a magnetic field profile similar 
to the density profile. Thereby t he predicted rotatio n measure (RM) profile agrees 
with the observed one (see Fig.l8)). lDolag et all (|200lh found a quasi linear correlation 



between two observables, namely the X-ray surface brightness and the RM r.m.s., both 
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Fig. 8 Comparison of RMs in simulations with observations of Abcll clusters, as a function 
of the distance from the cluster centre. The smooth lines represent the median values of 
|RM| produced by simulated clusters with masses above 5 X 10^"* Mq (upper line) and 3 X 
10^* M0 (lower line). The broken line represents th e median of combined data taken from the 
independent samples presented in lKim et al.l (Il991h and Clarke ct al. ( 2 00ll). We also in clude 
data (diamonds) for the three elongated sources observe d in A119 llFeretti et al.l I1999D , and 
for the elongated source observed in the Coma cluster (iFeretti et al.lll995l) . 



in observations and in tiieir MHD simula tions. By ex t ending the G ADGET2 code to 
follow the full set of ideal MHD equations. Ibolag et all (|2004l . l2005al ) performed several 
realizations of a cosmological volume which confirm these results at even much higher 
resolution. 

Recently, IBriiggen et al ] l|2005l ) confirmed all the previous results, which are based 
on SPH codes, with a simulation of the formation of a single galaxy cluster in a cosmo- 
logical framework, using a passive MHD solver implemented into FLASH, an adaptive 
mesh refinement (AMR) code (see Fig. |9]). 

Another interesting quantity to look at is the slope a of the magnetic field power 
spectrum (oc with fc being the wave vector). Within galax y clusters, a is predicted 
by the SPH simulations (|Dolag et al.ll2002l : [Ro"rdorf et al.|[2004l ) to be shghtly lower, but 
still very close to 11/3, which is the expected value fo r a Kolmogorov-like spectrum 
in 3D. The AMR simulation bv IBriiggen et al] ( 20051 ) nearly perfectly matches the 
Kolmogorov slope. 



4.4 Approximative simulations 

The fully consistent cosmological MHD simulations described above are very expensive 
in their requirements of computing resources; thus, for the investigat ion of various 
model s, simplified approaches are often used in astrophysics. Recently, iBertone et al.l 
( 20061 ) used semi-analytic simulations of magnetised galactic winds coupled to high- 
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Fig. 9 Slice through the centre of the simulated box hosting a massive galaxy cluster. The 
simulation used the AMR code FLASH and followed the evolution of a weak magnetic field 
seed. The fig ure shows the logarith m of the magnetic field strength, in Gauss, measured at 
z = 0.5. From lBruegen et al.l ll2005t) . 



resolution N-body simulations of structure formation to estimate limits for the fraction 
of the ICM which can be significantly magnetised. Interestingly, the fraction of the 
volume of the ICM, which is predicted to be magnetised at a significant level by the 
galactic outflows during structure formation, is very high. In fact, the volume filling is 
high enough to be compatible with observations and a moderately efficient turbulent 
dynamo, operating within the ICM, will easily amplify these fields up to the observed 
values. 

Additional to any dynamo action, further support for strong magnetic field ampli- 
fication during the process of structure for mation comes from the application of the 
Zel'dovich approximati on dZel'dovichl 1970l) to follow the MH D equations during the 
gravitational collapse i Bruni et al.l 200^ King fc Cole^l2006h . These papers demon- 



strate the existence of a super-adiabatic amplification due to the anisotropy of the 
collapse of the LSS within the cold dark matter paradigm. 

A novel aspect of including magnet ic field pressure into LSS simulations - even if 
in a simplified way - is investigated bv lCazzola et al" I l|2007l ). which followed a purely 
empirical approach. Inspired by the scaling between magnetic field and gas density, 
found in both simulations and observations i Dolag et al.ll200l[ ). the magnetic field is 
parameterised as B{z) = Bq{1 + zf"^ . Then an isotropic pressure (oc B^) term is 
added into the equation of motion. This approach allows a quite effective way to identify 
which models of magnetic fields, along with their parameters, provide a non-thermal 
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pressure support which is able to modify significantly the formation of structures. 
Thereby this approach gives a unique insight for the possible consequences of certain 
magnetic field scenarios. 



5 Cosmic Rays 

5.1 Modelling of Cosmic Rays (CRs) 

CRs may originate in a number of cosmological processes including jets due to mass 
accretion, large-scale shocks and supernova remnants (see Sect. [2]). One can conven- 
tionally distinguish two important classes of CR modelling in cosmological simulations. 

The first one is related to the direct modelling of the origin, spatial distribution and 
spectral evolution of energetic electrons and positrons, which are responsible for most of 
the non-thermal emission observed so far in cosmological objects. In the vast majority 
of the non-thermal lepton (e^) simulations, e^ are considered as test particles scattered 
and accelerated by electro-magnetic fields associated with the cosmological flows. In 
the highly relativistic energy regime, the leptons radiate very efficiently and therefore 
suffer from radiative losses. Most of the models are based on the linear kinetic equations 
for the electron distribution function without accounting for the back-reaction that the 
accelerated electrons have on the bulk plasma motions and on the MHD-turbulence. 

The second class of the CR simulations is related to the acceleration and propaga- 
tion of energetic baryons. The radiation inferred from the energetic trans-relativistic 
baryons is basically caused by their nuclear interactions with the ambient matter which 
produces 7-rays, neutrino and e^-pairs. Non-relativistic super-thermal baryons can 
produce some emission by inverse-Bremsstrahlung mechanisms. The most important 
feature of that class of processes, however, is that energetic non-thermal baryons could 
absorb a substantial fraction of the ram pressure of both non-relativistic and relativis- 
tic plasma flows. The back-reaction effects of energetic non-thermal particles on the 
plasma flows must then be taken into account. The most popular example of such a 
non-linear process is the d i ffusive particle accelera tion by astrophysical shocks (see e.g. 
iBlandford fc Eichleij[l987l : iMalkov fc DrurvlboOll . for a review). The process of diffu- 
sive particle acceleration by MHD shocks could be very efficient in highly turbulent 
cosmic plasma. 

In Fig.[l]of Sect.[2]we illustrated the simulated efficiency of the conversion into CRs 
of the power of shocks and of the bulk plasma motions in the particular cosmological 
context which represents the typical environment within galaxy clusters. The non- 
linear effect of the back-reaction of accelerated particles on large-scale plasma flows 
results in the specific temporal evolution of the particle spectra shown in Fig. [5] of 
Sect. [21 Moreover, as it was mentioned above, under a certain condition, thanks to 
the conversion of a fraction of the shock ram pressure into magnetic field energy, an 
efficient acceleration of baryons by MHD shocks in a turbulent cosmic plasma results 
in a strong amplification of the magnetic field in the shock upstream. 

The acceleration of relativistic electrons and their emission were modelled in dif- 
ferent cosmological contexts. To study the gradients of the spectral indices inferred 
from the r adio observations at 1.4 — 5.0 GHz of five luminous 3C radio galaxies, 
iMvers fc~Sp anglcr ( 1985, ) developed a model in which the radiation is due to an isotropic 
ensemble of rela tivistic electrons which ar e subject to synchrotron radiation losses (syn- 
chrotron aging) . iMvers fc Spangle ^ (|l985l) inferred the age of the electrons, namely the 
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Fig. 10 Results from a simulation intended to reproduce the merging system observed in 
A 3667. The left panel shows the simulated X-ray surface brightness (contours) and the radio 
emission at 1.4 GHz (grey scale). The image is 3.15 X 3.85 Mpc. The dashed and dotted lines 
refer to the location of the r adio spectral index (04 9) profiles displayed in the right panel. 
From lRoettiger et al.1 ||l999i). 



time since their acceleration, at various locations in the lobes, and consequently derived 
the speed of the separation of the hot spots from the lobe material. The inferred speeds 
were found to be in the ranee 10 000 - 30 000 km s" ^ These resuhs are consistent with 
the beam model, in which the lobe material is left behind by a hot spot advancing 
through the intergalactic medium at speeds of about 10 000 km s~^. 

lEnfjlin fc Gopal-Krishnal ( 200 ll ) modelled different scenarios for the evolution of 
the radio plasma inside the cocoons of radio galaxies, after the activity of the central 
engine has ceased. These authors discussed an analytical model for the evolution of 
a relativistic electron population under synchrotron, inverse Compton and adiabatic 
energy losses or gains. It was demonstrated that fossil radio plasma with an age of up to 
2 Gyr can be revived by compression in a shock wave of large-scale structure formation, 
caused during the merger events or the mass accretion of galaxy clusters. The scenario 
was applied to explain the origin of the diffuse radio emission found in clusters of 
galaxies, without any likely parent radio galaxy seen nearby. The model predicts the 
existence of a population of diffuse radio sources, emitting at very low frequencies with 
an ultra-steep spectrum, which are located inside and possibly outside galaxy clusters. 
Thereby, extended radio emission could be interpreted as tracing the revival of aged 
fossil radio plasma which is caused by the shock waves associated with large-scale 
structure formation. 

Some processes related to CRs were implemented in cosmological simulation codes. 
For instance, COSM OCR is a num erical code for the investigation of CRs in com- 
putational cosmology ( Miniatill200lh . The code includes a number of prescriptions to 
account for the diffusive shock acceleration, the mechanical and radiative energy losses 
and the spatial transport of the energetic particles into the cosmic environment. Pri- 
mary CR electrons and ions are injected at the shock sites according to the phenomeno- 
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Fig. 11 Results from a simulation of a shock interacting with a fossil radio plasma. The upper 
(lower) panels show t he time evolution of t he ga s density (magnetic field energy density) during 
the interaction. From lEntilin fc BriiggenI (|2002f ). 



logical thermal leakage prescription. Secondary electrons are continuously injected as 
a results of proton-proton (p-p) inelastic collisions of primary CR ions and thermal 
background nuclei. The code consists of a conservative, finite volume method with a 
power-law sub-grid model in momentum space. Note that in this numerical approach, 
the back-reaction of the non-thermal components (CRs and magnetic fields) caused by 
their pressure contribution to the thermal gas is neglected. 

To study the impact of CRs on galaxy and cosmic structure formation and evolu- 
tion, lEnfilin et all (|2007l ) developed an approximative framework which treats dynami- 
cal and radiative effects of CRs in cosmological simulations. These authors approximate 
the CR spectrum of each fluid element by a single power-law, with spatially and tempo- 
rally varying normalisation, low-energy cut-off, and spectral index. In this framework 
they included some approximate prescriptions for CR injection and acceleration by 
shocks, as well as CR transport and energy losses due to Coulomb interactions, ion- 
isation losses, Bremsstrahlung losses, and hadronic interactions with the background 
matter. These prescriptions are suited to be included into the global schemes of the 
numerical simulations of galaxy and structure formation. Although in such implementa- 
tion the description of the CR population is more simplistic than in the work described 
earlier, the dynamical influence of the CRs onto the underlying hydrodynamics is no 
longer neglected. This is not only important for the dynamics of the ICM itself but 
also for the injection of the CRs by shocks, which are altered by the presence of the 
non-thermal pressure support of the CRs. 



5.2 The quest for radio relics 

iRoettiger et all (Il999al ) were able to reproduce the main features of the extended pe- 
ripheral radio emission (the so-called radio relics) observed in A 3667 by combining the 
single merger simulations (Sect. l47l|) with a model for the in situ re-acceleration of the 
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Fig. 12 Results from a simulation of a galaxy cluster undergoing a major merger. The fig- 
ure shows the projected 'potential' radio luminosities for 1.13 Gyr old radio plasma, where 
PB/Pga,B = 0.01. For comparison, the bolometric surface X-ray luminosity (contours) is given. 
The total bolometric X-ray luminosity of th e cluster is 2 X 10 *^^ ergs~^ and the emission- 
weighted temperature is 3keV. Adapted from lHoeft et all (l2004f) . 



relativistic particles, as described in the previous section. In their models, they injected, 
into the ICM, relativistic electrons with a power-law spectrum, where the power-law 
index 7 — 3/(r — 1) -I- 1 is related to the gas compression ratio r — P2/P1 at the 
shock. They also related the age of the radio plasma fa to the distance d — KWgfa, 
using a weak-field/high-diffusion limit k = 1. By having effective shock velocities 
Us ~ 700 — 1000 kms~^ and by aging the synchrotron spectrum with the formal- 
ism of lMvers fc Spanglerl I igSSl ). they were able to reproduce the observed distribution 
of the spectral index for a magnetic field of ~ 0.6/iG at the position of the radio relic 
(Fig. llOfl . Since such configurations seem to be quite common in galaxy clusters, it 
naturally raises the question, why not all the clusters show such peripheral radio emis- 
sion. One possible explanation is that such shock structures are relatively short lived 
compared to the merger event itself. Also, it is necessary to have the presence of a 
large-scale magnetic field. It could further be that only massive clusters can provide 
enough magnetic field and strong enough merger events to trigger such a peripheral 
emission. 



To overcome this problem. lEnfilin fc Briiggenl l|2002l ) proposed that such radio relics 



could be made by a pre-existing fossil radio plasma illu minated by the shock waves that 
the merger events originate. In their models, following lEnfilin fc Gopal-Krishnal (|200lh 
which take into account synchrotron, inverse Compton and adiabatic energy losses and 
gains, they evolved the electron spectrum for the tracer particles, which represents the 
fossil radio plasma. Their simulation, using the ZEUS code, follows the evolution of a 
sphere of tracer particles hit by a shock front (Fig. [TTJ. Such a configuration nicely 
reproduces the filamentary radio emission and toroidal structures observed in many 
cases. These simulations also predict that the magnetic fields are mostly aligned with 
the direction of the filaments, as suggested by observational data. 
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iHoeft et"ai1 (|2004l ) investigated this id ea in a more realistic modelling by using the 
SPH code GADGET (jSoringel et al.ll200lh to simulate a merging galaxy cluster within 
a cosmological environment. Such a simulation showed that the probability for a shock 
wave to flare the radio plasma is highly suppressed in the central regions of galaxy 
clusters, unlike the peripheral regions, where they found illuminated structures of size 
up to Mpc scales (Fig. I12p . The reason for this is that the radio plasma ages much 
faster in the cluster centre than in the outer regions. In fact, in the cluster centre, the 
pressure of the radio plasma is higher and its energy losses due to the higher magnetic 
field are larger; moreover, the compression ratio of the shock wave is much higher in 
the low-density peripheral regions than in the cluster centre. 

It is worth noting that a necessary condition to form such relics is that the initial 
state of the fossil radio plasma is characterised by a ratio of the magnetic to thermal 
pressure, Pe/^gas, which has to be as low as 1 per cent to a llow shocks t o rev ive « 
1 Gyr old radio ghosts. It is also important to mention that iHoeft et al. I (|2004 ) find 
a high probability of radio emission, outside the shocks, due to the drained gas flows 
induced by the merger events, which transport material from the outskirts towards the 
higher density regions. Therefore in some cases the adiabatic compression seems to be 
sufficient to revive the fossil radio plasma. 

Recently, the picture of diffuse radio emission of galaxy clusters has changed, as 
more and more detailed observations indicate a more complex pic ture than the two 
main phenomena of radio halos and relics (see iFerrari et al. I boOSl - Chapter 6, this 
volume). It can well be that both scenarios described above are at work in galaxy 
clusters. The in situ re-acceleration of the relativistic particles at shocks might be 
responsible for the large, external relics, whereas the compression of the fossil radio 
plasma might lead to the smaller relics which are more disturbed in shape and are 
usually found closer to the cluster centre. 



5.3 The quest for radio haloes 

The accretion shocks from the LSS formation are also the location of the acceleration 
of the CRs that will be accreted into the LSS, spec ially within galaxy clusters. Using 
the COSMOCR code (,Miniatl. 20011 ). iMiniati et al.1 ([ioOli ) followed primary ions and 
electrons (injected and accelerated by diffuse cosmic shocks) and secondary electrons 
and positrons (produced in p-p inelastic collisions of CR ions with thermal ICM nu- 
clei) within a cosmological simulation. Under the assumption that the magnetic field 
produced by the battery effect refiects a fair representation of the true distribution 
of the relative magnetic field strengths within the LSS, they were able to predict the 
central radio emission (radio halos), mainly produced by secondary electrons in a self- 
consistent treatment (Fig. I13|l . 

Within their modelling they also reproduced the peripheral radio emission (radio 
relics), mainly produced by primary CRs, where the resulting morphology, polarization 
and spectral index match the observed counterparts. However, one has to note that 
the extrapolation of these simulations (which are on group scale) to the observed data 
(which are on the scale of massive clusters) might not be straightforward (Fig. I14|) . 
It is further worth commenting that even with an electron to p roton injection rati o 
^e/p of 10~^ (which is derived from observations, see Table 5 in lMiniati et al.ll200ll ). 
the predicted luminosities for the primary emission are still significantly higher than 
for the secondary one (Fig. I14|) : this prediction does not seem to be confirmed by the 
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Fig. 13 Predicted radio properties of halos within a cosmological simulation of a cubic vol- 
ume with a side length of 50 Mpc/h, resolved with 256^ grid cells. The left panel shows the 
synchrotron power of these halos at 1.4 GHz from secondary electrons. The right panel shows 
the synchrotron power from primary electrons, mostly located in the peripheral regions. Note 
that the values for the luminosities due to t he primary electrons should be scaled with the 
electron to proton injection ratio R^/p- From lMiniati et al.l ll200lh . 



observations. Although these results based on a secondary mo del generally agree wit h 
previous work based on simulations of massive galaxy clusters (|Dolag fc EnfiliijboOOl '). 
they highlight the general problem that one faces when it is applied the secondary model 
for relativistic electrons. This model predicts radio haloes for every galaxy cluster, 
which is in strong contradiction with current observations, as summarised in Fig. 1141 
which also i ncludes the upper lim its obtained by a recent observation campaign at the 
GRMT (see lBrunetti et al.|[2007l '). 



6 Turbulence 

6.1 Simulations 

During the growth of the LSS, galaxy clusters continuously accrete other structures, 
most of the time smaller structures (like galaxies), but sometime also objects with 
similar mass (major mergers). All these structures continue to move with supersonic 
or transsonic velocities inside the cluster. Together with the diffuse accretion (which is 
generally anisotropic) and the generation of turbulence by hydrodynamic instabilities 
induced by these bulk motions, the gas in clusters of galaxies generally contains an 
amount of kinetic energy which is not negligible compared to the amount of thermal 
energy. 

It is worth mentioning that different simulation methods reach good agreement 
in predicting that t he ratio of bulk ki netic energy to thermal energy is up to 15% in 
galaxy clusters (see Frenk et al."l999'). Based on cosmological hydrodynamic simula- 
tions, .Inoga mov fc Sunvacv (2003) pointed out that the broadening of the emission 
lines in the X-ray band (e.g. the iron K line) due to these expected bulk motions is 
appreciably larger than the broadening due to thermal motions. Thus, instruments like 
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Fig. 14 Total power of rad io halos observed at 1.4 MHz vs. cluster temperature. We plot the 
data from lLiang et all | |200C1|). which w ere partially rc-observcd by Feretti (2007, in prepa ration) 
together with data from lGovoni et al.l (2001); Bacchi ct al. ( 200a) : rVenturi et al.l ||2003D . Some 
additional upper limits are collected thanks to the help of F. Govoni. Additional upper limits 
fr om [Br unetti et al. (2007) are shown. We applied a secondary hadronic model, as described 
in lDolae fc EnfilinI 11200011 . to calculate the radio emission from the simulated galaxy clusters. 
We also added the predicti ons for the em i ssion from primary (upper line) and secondary (lower 
line) electrons taken from lMiniati et al.l ||2001| ) within the temperature range covered by the 
simulation. Note that the values for the luminosities for primary electrons should be scaled 
with the electron to proton injection ratio i?c/p- 



the ones designed for Suzaku or future instruments like XEUS will be able to infer such 
bulk motions from the analysis of the line shapes. 

Even after subtraction of the cluster peculiar velocity, large-scale bulk motions 
will be the main contributor to the deformation of the line shapes. The imprint of 
turbulence induced by hydrodynamic instabilities along flows (e.g. shear flows) might 
instead be more subtle to infer from the line profiles. 

In recent hig h-resolution SPH sim ulations of galaxy clusters within a cosmologi- 
cal environment, iDolag et al.l (l2005b[) used a novel scheme to treat artificial viscosity 
within the GADGET2 code (|Springei2005l ) to quantify how intense the shear flow must 
be to drive fluid instabilities. Being a quite common process within cosmic structure 
formation, these fluid instabilities driven by shear flows strongly increase the level of 
turbulence within the ICM (Fig. [T^ . 

Although this small-scale turbulence contributes a measurable effect to the non- 
thermal pressure within the ICM, the effect on the line shapes is much smaller than the 
broadening caused by the large-scale bulk m otions. Applying this method to a set of 
simulated galaxy clusters, IVazza et al. I (|2006l ) showed that for relaxed clusters the tur- 
bulent energy content in the ICM scales with the thermal one. The fraction of turbulent 
to thermal energy turned out to be rather independent of the mass of the system. Inter- 
estingly, this is consistent with the semi-analytical predictions of lCassano fc Brunettil 
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Fig. 15 The gas velocity field in a slice through the central Mpc region of a simulated cluster 
after subtracting the global mean bulk velocity of the cluster. The left panel shows a run 
with the original SPH viscosity, the right panel for a low-viscosity scheme. The underlying 
colour maps represent the ratio of turbulent kinetic energy to total kinetic energ y content of 
the par ticles, inferred after subtracting the local mean velocity, as described in iDolag et al.l 
d2005bl') 



( 20051 ). who estimated the amount of turbulent energy within the ICM in terms of the 
PAV work done by the infalling substructures. 



6.2 Observational evidence 

The spectroscopic resolution of current X-ray instruments is not high enough to mea- 
sure the line shapes and infer the level of turbulence present within th e ICM. However, 
using a mosaic of XMM-Newton observations of the Coma cluster, ISchuecker et al.l 
( 20041 ) were able to produce spatially-resolved gas pressure maps which indicate the 
presence of a significant amount of turbulence. Performing a Fourier analysis of the 
data reveals the presence of a scale-invariant pressure fluctuation spectrum in the range 
between 40 and 90 kpc which is well described by a projected Kolmogorov turbulence 
spectrum; at least 10 percent of the tota l ICM pressure is in t urbulent form. 

Alternatively, in a more indirect wav. lEnfjlin fc Vogtl ( 20061 ) argue that the recently 
reported Kolmogorov-like magnetic turbulence spectrum, which is inferred from Fara- 
day rotation measureme nts in the cool core of clu sters, can be understood by kinetic 
energy injection (see also lSubramanian et aL|[2006l ). Such a dynamo model predicts the 
correct magnetic field strength in cool core clusters for reasonable (but yet not directly 
observable) values of the hydrodynamic turbulence velocity and characteristic length 
scales. This result indicates that the magnetic fields might directly reflect the presence 
of hydrodynamic turbulence. Such models are directly related to the turbulence in- 
duced by buoyant radio lobes from the central radio galaxy which rise into the cluster 
atmosphere. On the other hand, Faraday rotation within non-cooling flow clusters with 
multiple extended radio sources (which therefore probe the magnetic field structure at 
different radii) can give alternative constraints on the magnetic field power spectra 



23 




0.06 



0.08 



0.1 

K [I Ape] 



0.2 



Fig. 16 The dots with error bars are the observed projected power spectral densities as 
inferred from X-ray observations of the Coma cluster, after subt raction of the shot noise; the 
dashed lines are model predictions. From [Schuecker et al. I (Hool). 



fsee iMurgia et aLll2004l : rGovoni et al.ll2006l ) and thus on the underlying hydrodynamic 
turbulence present in the ICM. 



7 Concluding remarks 



It seems that, in the last years, a consistent picture of magnetic fields in clusters of 
galaxies has emerged from both numerical work and observations. Simulations of in- 
dividual processes like shear flows, shock/bubble interactions or turbulence/merging 
events predict a super-adiabatic amplification of magnetic fields. It is worth mentioning 
that this common result is obtained by using a variety of different codes, which are 
based on different numerical schemes. Further support for such super-adiabatic am- 
plifications comes from an alytical estimates of the anisotropic collapse which use the 
Zel'dovich approximation ( Zerdovichlll97ol ). When this amplification occurs in fully 
consistent cosmological simulations, various observational aspects are reproduced (e.g. 
Fig. [SI; moreover, the final strength of the magnetic field reaches a level sufficient to 
link models that predict magnetic field seed at high redshift with the magnetic fields 
observed in galaxy clusters today. It is important to note that all the simulations show 
that this effect increases when the resolution is improved, and therefore all the current 
values have to be taken as lower limits of the possible amplification. 

Despite this general agreement, there are significant differences among the predic- 
tions on the structure of the magnetic fields coming from different models of mag- 
netic field seeds. In parti cular, there are seve r al differences be tween the up-scaled cos- 
mological battery fields ( Miniati et al.l boOll : ISigl et al] |2004| ) and the magnetic field 
predicted by the hig h-resoluti on si mulations of galaxy clusters that use either AMR 
( Briiggen etahl 120051 ) or SPH IPoTa g et al. 1999, 2002, 2005a). In the SPH case, it is 
possible to follow the amplification of the field seeds within the turbulent ICM in more 
detail. 
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A good visual impression can be obtained by comparing the regions filled with the 
high magnetic fields shown in Fig.[6]and|9l It is clear that the high magnetic field regions 
for the battery fields are predicted to be much more extended; this result leads to a 
fiat profile around the forming structure, whereas for the turbulent amplified magnetic 
fields the clusters show a magnetic field distribution which is much more peaked. Part of 
this difference originates from the physical model, as the cosmological shocks are much 
stronger outside the clusters. Somewhat less clear is how big it is the influence of the 
different numerical resolutions of these simulations to these discrepancies. Usually, the 
model parameters for such simulations can only be calibrated using the magnetic field 
strength and structure within the high-density regions of galaxy clusters. Therefore it is 
crucial to perform detailed comparisons with all the available observations to validate 
the simulations. Note that extrapolating the predictions of the simulations into lower 
density regions, where no strong observational constraints exist, will further amplify 
the differences among the simulations in the predictions of the magnetic field structure. 

Moreover, one has to keep in mind that, depending on the ICM resistivity, the 
magnetic field could suffer a decay, which is so far neglected in all the simulations. 
Furthermore, a clear deficiency of the current simulations is that they do not include 
the creation of a magnetic field due to all the feedback processes happening within the 
LSS (like radio bubbles inflated by AGN, galactic winds, etc.): this might alter the 
magnetic fleld prediction if their contribution turns out to be signiflcant. Also all of 
the simulations so far neglect radiative losses; if included, this would lead to a signif- 
icant increase of the density in the central region of clusters and thereby to a further 
magnetic field amplification in these regions. Finally, there is an increasing number of 
arguments suggesting that instabilities and turbulence on very small scales can amplify 
the magnetic fields to the observed nG level over a relatively short timescale. However, 
it is still unclear how such fields, which are tangled on very small scales, can be further 
processed to be aligned on large scales (up to hundreds of kpc), as observed. 

Concerning radio haloes and relics, the picture is only parti ally consi s tent (for a 
more detailed discussion of primary and secon dary models see iBrunettil ( 20041 ') and 
references therein as well as lFerrari et ahlboosl - Chapter 6, this volume). Radio relics 
seem to be most likely related to strong shocks produced by major merger events 
and therefore produced by direct re-acceleration of CRs, the so-called primary models. 
Although some of the observed features, like morphology, polarisation and position 
with respect to the cluster centre, can be reasonably well reproduced, there might be 
still some puzzles to be solved. On the one hand, direct acceleration of CRs in shocks 
seems to overestimate the abundance, and maybe the luminosity, of radio relics; on the 
other hand, simulations which illuminate fossil radio plasma can produce reasonable 
relics only starting from a small range of parameter settings. 

A similar situation arises for modelling the central radio emission of galaxy clusters 
by secondary models. The total luminosity seems to be reproduced using reasonable 
assumptions and also the observed steep correlation between cluster temperature/mass 
and radio power seems to be reproduced quite well. However, these models suffer from 
two drawbacks. The first one is that, in the framework of such models, every massive 
cluster produces a powerful radio halo, at odds with observations (Fig. I14|l . The other 
problem of secondary models is that the detailed radio properties are not reproduced. 
In fact, firstly, in most cases, the profile of radio emission is too steep, so that these 
models can almost never reprodu ce the size of the observ ed radio halos. Secondly, the 
observed spectral steepening (e.g. ICiacintucci et ahlboOSl ) cannot be reproduced. 
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One possible way out, at least for the first concern, would be that magnetic fields 
are extremely dynamical and transient and therefore light up the halo during a merger. 
However, there is also no indication from observations that clusters which show radio 
emission contain magnetic fields which are more intense than those in clusters without 
observable extended diffuse radio emission. On the contrary, the cluster A 2f42 has a 
magnetic field strength similar to the Coma cluster, but the upper limit on its radio 
emission is at least two orders of magnitude below the value expected from the corre- 
lation (Fig. I14|l . Note that both clusters are merging systems which are characterized 
by the presence of two central cD galaxies. This indicates that there should be further 
processes involved or additional conditions required to produce radio emission. It is 
worth noting that recent models, based on turbulent acce l eration, seem to overcome 
this problem (see IKuo et~aLl boosi: brunetti fc Blasil boosi : ICassano fc Brunettil boosl . 



and references therein). Cassano fc Brunetti (|2005h predicted that the probability for 
a galaxy cluster to show a giant radio halo is an increasing function of the cluster mass; 
these authors also reproduced the observed probability of ~ 30 % for a massive galaxy 
clusters to show extended radio emission. 

Being dominated by large-scale bulk motions, induced by the cosmological structure 
formation process, different simulation methods reach good agreement in predicting 
that the ratio of the bulk kinetic energy to thermal energy has an upper limit of 15 % 
in galaxy clusters (see lFrenk et al.ll 19991 ). However, the exact role of turbulence induced 
by hydrodynamic instabili ties within the com plex fiow patterns of galaxy clusters is 
still hard to quantify (see iDolag et aLlbOOSbl ). It should also be mentioned that the 
real viscosity of the ICM is hardly known, and can play a noticeable role within the 
ICM, as demonstrated in simulations w hich attempt to include p hysically motivated 
viscosity into a cosmological context (see lSiiacki fc Springellbood ). 

The dynamical role of CRs in galaxy clusters is not very well understood yet. Sim- 
ulations using different codes predict quite different relative pressure contained in CRs 
within galaxy clusters. Moreover, the simulations indicate that the relative importance 
of CRs in galaxy clusters also strongly depends on o ther non-thermal proc esses, like 
radiative losses and feedback from star formation (see IPfrommer et ahlbOOTl ). 

In general it has to be pointed out that, although the complexity of the physical 
processes treated by cosmological simulations have improved dramatically in recent 
years, this research field is still young. Therefore future work is expected to enlighten 
such complex processes within the ICM and will help to improve our understanding of 
the non-thermal components in galaxy clusters. Such work will also be needed to inter- 
pret the observational information which is overwhelmingly increasing in quality and 
size and is expected to become available soon with the next generation of instruments, 
particularly at radio band, like ALMA, LOFAR, SKA, EVLA and others. 
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